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Abstract 

We investigate the extended quasi-particle states in the mixed state of d-wave superconductors on 
the basis of the Bogoliubov-de Gennes equation. We prove that the quasi-particle eigen-states can be 
classified in terms of new topological quantum numbers which are related to the topological nature of 
the non-trivial phases of the quasi-particles. Numerical results for the quasi-particle eigen-states reveal 
the crossover behavior from gapless to gapped states as the flux density B increases. In the strong field 
region quantum oscillations appear in the excitation energy of the quasi-particles. 

1 Introduction 

The quasi-particles in the mixed state of d-wave superconductors has attracted much attention in connection 
with the peculiar low-energy excitations under a magnetic field in high- Tc cuprates [|l|, ||, |[ ||, ||, 0, ^, 



^ 10, |11|. In this paper we investigate the extended quasi-particle states in extremely type-II d-wave 
superconductors on the basis of the Bogoliubov-de Gennes (BdG) equation, focusing on the topological 
nature of the quasi-particles. It is well-known that the order parameter has a non-trivial phase in the 
presence of vortices, that is, the phase (/)(r) is a multi-valued function satisfying the relation, V x V0(r) = 
2TTj2i S{r — 'Ri)ez, in the mixed state where vortices (|| e^) are located at Ri, R2, • • • , Rat. Then, the 
BdG equation contains a multi-valued function, which makes the equation intractable. Anderson noticed 
that the phase (/)(r) in the BdG equation can be eliminated in terms of a singular phase transformation, by 
which the BdG equation is transformed into the one containing V(j){r) Q. This is an important observation, 
because V0(r) is a single- valued function, that is, the transformed BdG equation includes only single- valued 
functions. However, it turns out, as shown in this paper, that the topological character of the system is not 
fully incorporated into the wave functions by this transformation. One should note that a superconductor 
containing vortices may be considered as a multi-connected system. Then, the quasi-particles moving around 
vortices show the AB effect [|ll[ . From this fact one understands that the wave functions of the quasi-particles 
should be path- dependent single-valued functions in the presence of vortices. 

In this paper we clarify the topological character of the non-trivial phases of the wave functions satisfying 
the BdG equation and show that the quasi-particle eigen-states have a "hidden" quantum number which 
is related to the homotopy class of the classical orbits of the quasi-particles. We also develop a theoretical 
scheme for solving approximately the BdG equation in the field region Hd <C <C Hc2, correcting the 
incomplete treatments for the topological singularity of (/)(r) in previous works. Numerical solutions for the 
extended quasi-particle states in d-wave superconductors are briefly presented. It is shown that the gapless 
quasi-particle state changes to a gapped one with increasing the flux density. In the strong field region 
quantum oscillations appear in the excitation energy. 
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2 Formulation 



Consider the BdG equation for a d^y-waYe superconductor, 



/io(r,V) A(r,V) 
A*(r,V) ~Kir,V) 



Ua(r) 



E 



Va (r) 



where 



and A(r, V) is the gap operator defined as 



A(r,V) 



fcl 



■{a.,{a„A(r)}}-^A(r){a,,9,}0(r), 



(1) 
(2) 

(3) 



with {A, B} = ^{AB + BA) In eq.(3) 0(r) is the phase of the gap function A(r) and the second term 
is added for recovering the U(l)-syniinetry |l3j. In this paper we investigate the extended quasi-particle 
states outside the vortex cores in the field region Hd <C ^ i?c2, so that we may approximate the gap 
function as 

A(r)~ Aoei^(^), (4) 

neglecting the spatial dependence of the gap amplitude. In the case where N vortices (|| z-axis) are located 
at Ri,---,Rn in the xy plane, the phase <j) in eq.(3) is dependent on the vortex positions, i.e., 4>{v) = 
(^(r; Ri, • • • , Rn), and is a multi- valued function satisfying the relation. 



N 



V X V0(r; Ri, • • • , Rn) = 27r ^ 5{v - Ri)e^ 



(5) 



i=l 



To express the BdG equation in terms of only single- valued functions, removing the multi- valued function 
^(r), Anderson proposed a transformation for the wave functions Ua and in eq.(l) as follows H, 



(6) 



This transformation has been utilized by several authors for studying the low energy excitations in the mixed 
state of d-wave superconductors 0. On the other hand, Franz and Tesanovic pointed out that the phase 
can be eliminated also by the transformation, Ua e^'^^u^, Va e~^'^'^Va, ii (j>e + fph = <t> l^j 0' ^^'^ phase 
factor e^'^ in eq.(l) can really be excluded from the BdG equation by these transformations, but we claim 
that these manipulations should not be considered as the transformations. Let us suppose that the wave 
functions in the A^- vortex state, MQ,(r), ^^(r), are expressed as 



ua(r) = 'S„^(r)ei(3-/^)'^('-;R-i.--.R-N) 
fa(r) = {;„^(r)e-i(H/^)'#'(r;R-i.-,R.w) 



(7) 



when the non-trivial phases in the wave functions are explicitly extracted. In eq.(7) /i is assumed to be a half 
and integer^ namely, /i = ±i,±|,---. Note that the phases in eq.(7) satisfy the Franz- Tesanovic condition 
mentioned above ((/)e + (ph ^ <!>) and eq.(7) is also reduced to the Anderson's transformation (6) for /i = ±i 
if eq.(7) is considered as the transformation, {ua,Va) — * {uonVa)- Hence, it is clear that the phase factor e'^ 
is eliminated in the equation for UonVa- The new equation for Ua and Va contains explicitly a number /i, so 
that the energy eigenvalues explicitly depend on /i. Furthermore, it turns out from the calculations shown 
later that the energy eigen-states can be classified in terms of /i. From these results one can interpret /i as 
a quantum number specifying the quasi-particle eigen-states in the A^-vortex state. Then, a set of quantum 
numbers a in eq.(7) is expressed as a = (n, fj,), where n denotes symbolically the other quantum numbers. 

To be convinced of our statement, let us now study the topological character of the wave functions. We 
introduce the operator 7'(C) which moves a function along a path C jl^. Let us take a closed path Ci 
around the ith vortex at R^ for C. Suppose that Ci does not include any other vortices inside it (see Fig.l). 
Then, the operation of V{Ci) on the phase (/)(r) induces the shift. 



P(C,)0(r : Ri 



(r :Ri,---,RN)-f 27r. 



(8) 
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Suppose that the wave functions, Ua{r) and Wa(r), have non-trivial phases, Qu{r) and Qy{r), i.e., Ua{r) = 
{ia(r)e"'®°'-''' and Vaiv) = Va{r)e^^'^'^''\ Since ^^(r) and Va{r) should be single-valued functions, we have in 
general 

V{Q)uo^=uJ^-', V{C,)v^=vJ^^"', (9) 
with £ and m being integers. Then, eq.(5) indicates the relations. 



(10) 



under the conditions, V{Ci)ua — Ua and V{Ci)va = Va- Furthermore, from the gap equation, A(r) 
faUaV^e^^^''~'^°'^ (scc Appendix A), one may assume the equality, 



Then, from the operation of V{Ci) on both sides of eq.(ll) it follows. 
Comparing this relation with eq.(8), one finds 



1 



m = l, or £= — — 11. m = — f- + /u), 
with being a half and integer. Since these relations are simultaneously satisfied when 



e^(r)= (i-M)0(r) 
0a(r)= -(i+M)0(r) 



(11) 
(12) 

(13) 
(14) 



we obtain eq.(7). We can also show that the gap equation expressed in terms of the wave functions (7) yields 
4>(r) for the phase of the order parameter (see Appendix A). Note that eq.(7) holds for = 1, i.e., the single 
vortex state. As is well-known, the wave functions in the 2D single vortex state are expressed as 



u„^(r,0) = {i„^(r)ei(^-'^)« 



(15) 



where ir,9) is the polar coordinates jl^, which indicates that the phase 4'{r) is expressed in terms of the 
geometrical angle 6, i.e., 0(r) = 9, in this state and then n (= ±i, |, • • •) accords with the quantum number 
for the 2D angular momentum, though in general /i is a topological one. Thus, our conjecture (7) is consistent 
with the well-known result in the single vortex state. Note that /i is still a rigorous quantum number in the 
system with a curved vortex, in which the angular momentum is not conserved. 

Let us now study the BdG equation (1), using eq.(7). In the present d-wave case we must be careful 
about the order of differentiations in eq.(3), because dxdy(p{r) ^ dydx4>(j) in the vortex state To 
avoid the ambiguity in the order of differentiations we symmetrize the differential operations in eq.(3) as 
SxSye*'^^"") {a:r,9y}e*'^(""). Then, substitution of eq.(7) into eq.(l) yields 



/io(r, V;a*) 
n*(r,V;-/.) 

where 

/io(i", V;^) 
Vs is the superfluid velocity defined by 





'^71 fJ, 


= E 






f^7l fJ, 







-I 



m 



m tic 



and the off-diagonal component, n(r, V; /x), is given as 

Ao 



n(r,V;M) 



• dy(j) - i^l{^x,^y}(t) - \^i{dy(j) ■ dx + dx(j) ■ dy) + d^dy 



(16) 



(17) 



(18) 



(19) 
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To perform explicit calculations we introduce an approximation in the following. We investigate the region 
Hci <C -ff ^ Hc2 in the extremely-type II superconductors of k 3> 1, i.e., ^ ^ A^. In this field region one 
can assume the relation, ^ <g; d <^ Xl, where d is the lattice constant of the flux-line-lattice. The spatial- 
dependent flux density in this region, B(r), can be split into the average value, Bav, and the deviation from 
it, <5B(r), i.e., B(r) = Bav + ^B(r). Since the lattice constant d is much shorter than the London penetration 
depth Ai in this field region, the vortex currents flowing around vortices are heavily overlapped and, as a 
result, canceled out, which indicates the relation, |(5B(r)| ^ |Bav|- Consequently, the spatial-dependent term 
(5B(r) is neglected in the 0th order approximation, B(r) ^ Bav In this paper we confine ourselves to the 
lowest order approximation. Then, in this approximation we have the relation from the Maxwell equation 
as follows, 

Aire* 

VxB(r) = p,(r)v,(r) ~ V X Bav = 0, (20) 

c 

where Ps(r) is the local superfluid density. Since /Os(r) ~ const. 7^ outside the vortex cores, one can assume 
Vg ~ in the present approximation. Then, noting eq.(18), we find the approximate relation, 



V0(r) H£A(r) ^ |^Ao(r), 
he he 



(21) 



where the vector potential Ao(r) is related to the average flux density as Bav = V x Ao(r). Since Ao(r) = 
(0,i?a;,0) with B = |Bav| in the Landau gauge {B < H), we arrive at the explicit expression for the 
non-trivial phase as 

2e 27r 

d,<j> = 0, dycl) = —Bx = —Bx, (22) 

he (po 

where (f>o is the unit flux, (po — Eq.(22) indicates the relation. 



dxdy(l) = ^B, 
he 



dydxii = 0, 



(23) 



that is, the differential operations for do not commute, dxdycj) ^ dydx4>, as expected. Then, we have the 
topological singularity in the present approximation, 



(V X Vcj){Y% = {dxdy - dydx)(t>ir) = 2Trm 



(24) 



noting B — m(j)Q with m being an integer. This result corresponds to that in the "continuity-approximation" 
for the singularity, i.e., (V x V^(r))^ ~ 27r(X;f S{r - Ri))av. = ^nm 

Let us now solve the BdG equation, using the above approximation. Under the approximation eq.(16) 
can be rewritten in the 2D case as follows, 



ho{P,Q) 
n*(P,Q;-^) 



U{P,Q;^l) 
-ho{P,Q) 







= E 













where 



and 



n(P,Q;/^) = ^nLo,[i^i + 2^PQ\. 
Zef h 



Here, = ^ is the cyclotron frequency and 

Q={^—By~^Px, 

^ e ' 



P 



u — Bx) , 

c ' 



(25) 

(26) 
(27) 

(28) 



withpa = —ihda- Note that Q and P satisfy the commutation relation, [Q, P] = ih. Thus, the BdG equation 
is expressed in terms of only Hermitian operators and is gauge-invariant, we remark that symmetrizing the 
differential operations in eq.(3) is essential for obtaining the above physically-acceptable results. Now let us 

introduce the raising and lowering operators, and a, as Q = ^^^(^ " '^^) ^^'^ ^ ^ \/§('^ + '^^)' noting 
the above commutation relation for Q and P. Then, eq.(25) is rewritten as 



2/i(a^a 

2eF 



^/i(a' 



-2/i(ata + 



















TlLOc 





(29) 
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where v'^^ — —ivnfi and e'p = To solve eq.(29) numerically it is convenient to introduce the expansion 
in terms of the Landau levels, 

oo oo 

Unt, = a'afe+l I ' ^'np = "2fc+2 I ^) ' (30) 
fc=0 k=0 

where |fc) is the fcth Landau level, i.e., a\k) = ^/k\k — 1) and a^\k) = ^/k + l\k + 1). It is seen that the 
off-diagonal components in eq.(29) induce the mixing between the Landau levels, that is, the Landau levels 
are not the eigen-states in the present d-wave case. Eq.(29) also indicates that the quasi-particle eigen- 
states in the mixed state of 2D d-wave superconductors can be classified in terms of two quantum numbers, 
namely n, which has the topological origin, and n, which is reduced to the Landau-level index in the normal 
state (Ao — > 0). The quantum number /i is a rigorous quantum number in the mixed state of type II 
superconductors, which is independent of approximations used in this paper. From the above results we 
finally obtain the coupled equations for the coefhcicnts, , as 

\^,AoVik-mk-l)aZ, + [Mfc - e'^K'^ - l^.Ao^ (k + l){k + 3)aZ, = ^^a^. (31) 
for /c = 1, 3, 5, • • •, and 

- i/iAov/(fc-4)(fc-2K^:5 - [M(fc - 1) - e'^WZ + l^,A„^^k(k+^H'+3 = (32) 
for /c = 2,4,6,- • •, with Aq = jf. 

3 Numerical results and discussions 

Let us now present numerical solutions of eqs.(31) and (32) obtained by a numerical diagonalization method. 
Fig. 2 shows the positive energy eigen- values for the quantum numbers, /x = 5, |, §• The parameter values in 
these calculations are chosen as e^^ = 200 meV, Ag = 10 meV, and the free electron mass We is assumed in 
the cyclotron frequency, i.e., Wc = '^rH' seen in this figure, we have two quasi-particle eigen-states with 
nearly-zero eigen-values. These states appears in the case of Tiuoc <C Aq in our calculations. The existence 
of the zero-energy state has been predicted in the approximate calculations in ref.[7]. These zero-energy 
states are formed mainly by the quasi-particles lying in the nodal directions in the fc-space. It is also seen 
that the two or four consecutive levels have close eigen-values in the low energy region. This result comes 
from that the state with index k does not couple with that with fc -I- ^ of ^ < 2 or £ < 4, as seen in eqs.(31) 
and (32). In Fig. 3 we plot the energy eigen-values in the A* = 5 band for three values of B. Note that the 
relation hujc Aq holds in these cases, since TilOc ~ 1.2meV^ Aq for _B = 10 T in the present choice of 
the parameter values. Then, the energy levels are essentially the same in this field region, except that the 
dispersion becomes steeper as B increases. 

To investigate the crossover behavior due to the variation of the flux density B from TllOc ^ Aq to 
hujc ^ Aq we calculate the lowest excitation energy as a function of B, choosing a small value for Aq. Fig. 4 
shows the i3-dependence of with n= 1, n — ^ in the case of Aq — 2 meV. In this choice of the parameter 
values we find HlOc = Aq for B ~ 16T. As seen in this figure, the lowest energy eigen-state has zero-energy in 
the weak field region of huJc ^ Aq, but it abruptly changes to a gapped state above some value of B, which 
indicates that the quasi-particle eigen-states acquire a full energy-gap in the strong field region, though the 
symmetry of the gap function is the same as in the Meissner state. This result may be intuitively understood 
in the following way. The quasi-particles lying in the nodal directions with Px ^ Pf cannot turn by the 
effect of a weak Lorentz force when the pair-potential barrier is high. Then, in the weak filed region the 
quasi-particle states in the nodal directions are expected to be gapless. However, in the strong field region 
where the Lorentz force is strong enough, i.e., huJc ~ Aq. the quasi-particles in the nodal directions can run 
through the pair-potential barrier. In order to excite the quasi-particles going through the barrier of the 
pair-potential we need a finite excitation energy, which leads to a finite energy gap. It has been proved in the 
GL limit, in which the Landau-level splitting is completely neglected, that the d^jy-component is induced in 
the mixed state of dx2^y2-wave superconductors [ p^ , which also creates fully-gapped quasi-particle states. 
However, the origin of the energy gap in the GL limit is completely different from that in the present theory. 
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This result seems consistent with the experimental observation for the thermal conductivity in the mixed 
state Q. In the strong field region of huJc > Aq, namely B > 16T in the present case, the lowest excitation 
energy shows oscillatory behavior in the -B-dependence, as seen in Fig. 4. These oscillations reflect the fleld 
dependence of the diagonal components of eq.(29), that is, they are the quantum oscillations which lead 
to the de Haas-van Alphen effect. Thus, our theory for the extended quasi-particles in the mixed state of 
d-wave superconductors can describes the crossover behavior from the gapless phase in the weak field region 
to the gapped phase showing the quantum oscillations in the strong field region. 

In summary we investigated the extended quasi-particle states of d-wave superconductors in the field 
region, Hd H <^ Hc2, on the basis of the BdG equation. We found new topological quantum numbers 
which classify the quasi-particle eigen-states in the mixed state of type-II superconductors. The phase of 
a qausi-particle should be path-dependent in the presence of vortices, which is essentially the same as in 
a multi-connected system showing the Aharnnov-Bohm effect. Hence, one understands that the quantum 
numbers carrying information about the homotopy class of the orbit of the quasi-particle appear in the 
mixed state of type II superconductors. We also presented numerical solutions of the BdG equation in the 
field region, Hd H Hc2- It was shown that the quasi-particle states show the crossover behavior from 
gappless to gapped states as the flux density increases. In the strong field region of hiVc > Aq the quantum 
oscillations appear in the excitation energy. 
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A Appendix 

In this Appendix we derive the gap equation in the present daj^-superconductors. The gap function in 
unconventional superconductors has generally the form, 

A(r,r') = y(r-r')(V'T(r)V'i(r')) 

= y(r - r') J2{MrK{r') [l - ME^)] - v* ir)u.{r') ME,)] , (33) 

■i 

where V(r — r') denotes the interaction between superconducting electrons. We introduce the Fourier trans- 
formation for the relative coordinates, s = r — r' as 

A(r, r') = A(R, s) = ^ A(R, k)e'^-^, A(R, k) = / dsA(R, s)e-i''■^ (34) 

k 

with R = i(r + r'). The wave functions, Ui{r) and Wi(r), in eq.(33) are expanded in powers of the relative 
coordinates as 

M,(R + is) = w,(R) + i(s • Vr)u,(R) + i(s • Vr)2u,(R) + • • • , (35) 
<(R - is) = <(R) - i(s • Vr)<(R) + i(s • Vr)^<(R) + • • • , (36) 

For 

8- dxywave superconductor one can utilize the following functional form for the interaction, 

^(«) = § E = ^dfdfS{s). (37) 

Then, substituting eqs.(35), (36) and (37) into eqs.(33) and (34), we find the Fourier component of the gap 
function, 

A(R,k) = i|^ y'ds5(s).5f9f{l?,(R,s)e-''-^}[l-2/F(i?0], (38) 

where 

A(R,s) = lis ■ Vr)'m,(R) • <(R) + lu,(R) ■ (s • Vr)'<(R) 
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-(s.VR)ii,(R)-(s.VRK(R) 



(39) 



Assuming A(R, k) ~ Ao{Il)kxky in the present d-wave case, we extract the (ia;y-component from eq.(38) as 
follows, 

A(R,k) = -i^^[a^a^I?,(R,s)]^^Jl - 2ME,)]kJy 
Thus, the gap function is obtained as 

Ao(R) 

where 



(40) 



(41) 



[9:a^2?.(R, s)]^^^ = i{a«a>,(R) • <(r) + u,(r) • a«9>;(R) 
- u.(R) • a>*(R) + a,«7.,(R) • <(r)] }. 

Note that the above gap equation is not invariant under the U(l)-transformation, 

Ao(r) ^ Ao(r)e'^('-\ «„(r) ^ w„(r)e^'xW, v^{v) ^ v^{v)e-"-'^^'\ 



(42) 



(43) 



The U(l)-symmetry may be recovered by replacing the derivatives in (A. 10) with "covariant" derivatives, 
da ^ da — ^4>a wherc (f>a = da4>, 4> being the phase of Ao(r), i.e., Ao(r) — |Ao(r)|e"^^""^ [|l^. Furthermore, 
introducing the symmetrization for the differential operations, dxdy — + {9a;, dy}, we rewrite the gap equation 
(A. 9) as follows. 



^o(r) = ^ \{dx - ]^(t)x,dy - ^0y}u,(r) • v*{r) + u,{r) ■ {d^ - ^(j)^, dy - ^'?!>y}<(r) 



{dx - ]^(t)x)u,{Y) ■ {dy - ^0j,)w*(r) - [dy - ^0y)wj(r) • [d^ - ]^(t)x)v*{v) 



tanh ■ 



E, 



(44) 



2kBT 

Here, we change the notation, R — s- r. Suppose that the wave functions, Ui and Vi, are expressed in terms of 
the non-trivial phase ^(r) as 



u„^(r) = u„^(r)e'(^-f)*('-) 
^^„p(r)=«„^(r)e-'(HM)0(r) ' 



(45) 



with i — in, fi). The new wave functions, 'Un/i(r) and Vn^ij), are assumed to have only trivial phases. Then, 
substituting eq.(45) into eq.(46), we obtain the gap equation in terms of Unp, and Vn^,, 



Ao(r) = -e'^C-) • ^ Yl [i^^ " f ^« " J'f>y}<^^i^) ' 

F 



+M„^(r) • {dx + ^-^<t)x,dy + y0y}w;^(r) - (9^; - y'/':r)unp(r) • (9y + y (/'y)w;;^(r) 



{dy - ^-^(t)y)u^i,{Y) ■ {dx + y0x)w,V(r) 



tanh ■ 



which indicates that the gap function Ao(r) have the non-trivial phase, (/)(r). 



(46) 
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in which the terms of G ^ are neglected one gets the approximate result given in eq.(24). 

[19] T. Koyama and M. Tachiki: Phys. Rev. B 53 (1996) 2662. 
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Figure 1: Closed path d around the vortex located 
at Ri. Solid circles denote the vortex cores. 



Figure 3: Energy eigen-values of /i = i-band for the 
flux densities, i? = 1, 5, lOT. The relation TilOc <C Aq 
is satisfied in these flux densities. The zero-energy 
state appears in this case. 
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Figure 2: Energy eigen-values for the quantum num- 
bers = ^ , I , I at i? = 5T in the low energy region. 
The eigen- value, has zero- value within numeri- 
cal errors. 



Figure 4: Dependence of the lowest excitation energy 
on the flux density B. In the weak fleld region, where 
hujc Aq is satisfied, the zero-energy quasi-particle 
state appears. However, it changes to a gapped state 
above some value of B. In the strong field region 
satisfying the relation huc > Aq the quantum oscil- 
lations appear in the quasi-particle excitation energy. 
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